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Dependence of fluorescence-level statistics on bin time size in a few-atom 

magneto-optical trap 

Sungsam Kang, Seokchan Yoon, Youngwoon Choi, Jai-Hyung Lee, and Kyungwon Arjj 

School of Physics, Seoul National University, Seoul 151-742 Korea 

(Dated: Received February 2, 2008) 

We have analyzed the statistical distribution of the fluorescence signal levels in a magneto-optical 
trap containing a few atoms and observed that it strongly depends on the relative size of the bin time 
with respect to the trap decay time. We derived analytic expressions for the signal distributions 
in two limiting cases, long and short bin time limits, and found good agreement with numerical 
simulations performed regardless of the size of the bin time. We found an optimal size of the bin time 
for minimizing the probability of indeterminate atom numbers while providing accurate information 
on the instantaneous number of atoms in the trap. These theoretical results are compared with actual 
experimental data. We observed super-Poisson counting statistics for the fluorescence from trapped 
atoms, which might be attributed to uncorrelated motion of trapped atoms in the inhomogeneous 
magnetic field in the trap. 

PACS numbers: 32.80.Pj, 42.50.Ar, 02.60.Pn 



I. INTRODUCTION 

One of the long-sought experimental capabilities in 
modern atomic physics and quantum optics is the ability 
to load a single atom in a microscopic volume for an ex- 
tended time and to manipulate and probe its internal and 
external states at will. In recent years, several groups 
have developed technics for trapping and controlling a 
single or a few neutral atoms based on tightly localized 
magneto-optical traps or dipole traps p], 0, LJ 3 • Sin- 
gle or a-few-atom traps have been applied to wide range 
of fields such as cavity quantum electrodynamics studies 
0, experiments on single-photon generation on demand 
p , and even archeological dating of ancient aquifers Q . 

The most distinctive signature of single atom trapping 
is the quantized fluorescence signal. When the number of 
trapped atoms is decreased to single-atom level, the flu- 
orescence signal from atoms exhibits stepwise underlying 
variation in time and the size of the fluorescence signal 
with respect to a background level is interpreted as be- 
ing proportional to the instantaneous number of atoms 
in the trap. Such stepwise fluorescence signals have been 
regarded as the most definitive evidence for single atom 
trapping. With this understanding, one can obtain the 
atom-number distribution in the trap from the histogram 
of the fluorescence signal levels and can also identify indi- 
vidual loading and loss events of atoms in the trap |j, Q • 

In actual single- atom trap experiments, since the fluo- 
rescence signal from a single atom is extremely weak, one 
needs to choose a bin time for photon counting as long as 
possible in order to achieve enough signal to noise ratio. 
If the bin time is too long, however, the atom number 
can change several times during the bin time and thus 
the observed fluorescence no longer provides accurate in- 
formation on the instantaneous atom number. 
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From the experimental point of view, therefore, sev- 
eral questions naturally arise regarding the conditions 
under which the fluorescence measurement should be per- 
formed: what will be the optimal size for the bin time, 
what determines the shape of the signal distribution and 
thus what information one can get from the observed sig- 
nal distribution. The purpose of this paper is to answer 
these questions. 

This paper is organized in the following way. In Sec. 
II, we define the problem and derive analytic expressions 
for the signal distributions in two limiting cases, long and 
short bin time limits, along with signal-to-noise consider- 
ations. These results are compared with numerical simu- 
lations in Sec. Ill, where an iterative method and Monte 
Carlo simulations are employed to calculate steady-state 
atom number distribution functions and the signal distri- 
butions regardless of the size of the bin time. In Sec. IV, 
an optimal size of the bin time is identified for minimizing 
the probability of indeterminate atom numbers while pro- 
viding accurate information on the instantaneous num- 
ber of atoms in the trap. The analytic expressions and 
numerical results are then compared with experimental 
results in Sec. V. It is demonstrated that the experimen- 
tal signal distribution is well fit by our theoretical model 
and from observed signal distributions one can extract 
information not only on the number of atoms but also on 
the state of atoms in the trap. In Sec. VI, we summarize 
our findings and draw conclusions. 



II. THEORETICAL CONSIDERATION 

In a few- atom trap, the fluorescence signal from atoms, 
induced by a probe laser or by a weak trap laser itself in 
the case of a magneto-optic trap (MOT), is proportional 
to the number of atoms in the trap. The fluorescence sig- 
nal is measured with a photodetector, usually in photon 
counting mode with photon counting electronics. Sup- 
pose the signal counts are successively taken in time for 



a preset bin time of At. The signal counts Si measured 
in i th time bin, specified as t{ < t < ti+\ with ti = iAt 
(i = 0, 1, 2, . . .), can be written as 



Si 



f 



U+At 




(i) 



where N(t) is the instantaneous number of atoms in the 
trap (N(t) = 0, 1, 2, . . .), a,j(t) is the counting rate of flu- 
orescence from j th atom, and b(t) is the counting rate of 
background signal such as detector dark counts and scat- 
tered laser light or stray room light. The bin time At is 
assumed to be much larger than spontaneous emission 
lifetime of atoms, typically tens of nanoseconds. The sig- 
nal Si is truncated to the nearest integer by the counting 
electronics. 

The instantaneous number of atoms N(t) rapidly fluc- 
tuates due to various stochastic processes. Temporal 
change of its probability distribution function P/v(£) is 
governed by the following master equation: 



dP N 
dt 



RPn- 



R + T 1 N + T 2 



N 



~^N 



+i\ (TV + i)p N+1 + r 2 f N + 2 ] p N+2 , (2) 

where N = 0, 1,2, . . . with P_i = 0, R is loading rate 
of atoms into the trap, Ti is the one-atom loss rate due 
to collisions with background gas, T 2 is the two-atom 
loss rate due to light-assisted intra-trap collisions Q- 
The master equation Eq.(|2j) cannot be solved analyti- 
cally. However, for a microscopic trap with only a few 
atoms in a volume of a few micron in diameter, the two- 
atom loss terms proportional to T 2 are negligibly small, 
and thus approximate expressions for Pjsr and the num- 
ber correlation function (N(t)N(t + r)) t can be obtained 
|9| with ( ) t denoting a time average. For now, we just 
neglect the two- atom loss terms. The analysis including 
these terms will be discussed later. 

Without the T 2 terms, the master equation Eq.© be- 
comes the same as the simple birth-death model of popu- 
lation U3 , yielding a Poisson distribution in steady state 
with a mean value N = R/Y\ and a variance a 2 = N. 
For a few- atom trap with N ~ 1, the correlation decay 
time r of the number correlation function is given by 
r = 1/ri and is the measure of the average time dur- 
ing which N(t) remains constant [3|. When T 2 terms 
are not negligible, the correlation decay time is a compli- 
cate function of R, Ti and T 2 and always smaller than 
1/ri due to the additional two-atom loss process 9]. We 
denote the correlation decay time in this case as r e ff in 
order to distinguish it from the above definition of r for 
the T 2 = case. The correlation decay time is also called 
the trap decay time for macroscopic traps. 

The signal distribution shows much different behaviors 
depending on the size of At with respect to r. For an- 
alytic analysis we consider two limiting cases, long bin 
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FIG. 1: Number of atoms in the trap fluctuates in time due to 
stochastic processes. Two limiting cases of bin time, At ^> r 
and At <C r, for the integral in Eq. JJJ are shown. 



time (At ^> r) and short bin time (At <C r), as illus- 
trated in Fig. E] In the long bin time limit, there exist 
many loading and loss events in a single bin time whereas 
the atom number hardly changes during the bin time in 
the short bin time limit, and thus these two limiting cases 
lead to quite different signal distributions. 



A. At > r limit 

Assume that N ~ 1 <C ajAt, bAt and that dj(t) and 
b(t) fluctuate much faster than r. Under this assumption, 
for a given N the variations in ctjAt and b At go like 
square root of those, respectively, and thus much smaller 
than a,jAt and bAt themselves. On the other hand, the 
variation due to N change is as large as ajAt. Therefore, 
in evaluating the integral in Eq. JJJ, we can neglect the 
fluctuations in ctj and b and replace them with their mean 
values a and 6, respectively. 



ti+At 



[N(t)a + b]dt = (Nia + b) At , (3) 



where 



i pti+At 

Ni = — I N(t)dt 



At 



(4) 



u 



is the time-averaged atom number in the i time bin. 
Since At ^> r, N^ which is no longer an integer, fluctu- 
ates around the mean value TV with a new variance <j 2 , 
which is not the same as the variance a 2 of P/v distribu- 
tion above. In fact, a 2 — > as At — > oo. The probability 
distribution P(Ni) can be obtained by the central limit 
theorem as a Gaussian distribution, 



P(Ni] 



1 



27ra 



exp 



•i^ 



Nf 



(5) 



where the variance a is proportional to the original vari- 
ance a 2 , which equals TV for a Poisson distribution, and 



inversely proportional to the sample size, which is in the 
order of At/r. The exact calculation for a Poisson dis- 
tribution (T2 = case) is given below: 



a 2 



(Nf) - (N z y 

rti+At 



rU+At 



= ^| dtj t dt'N(t)N(t'))-N 2 

-1 nAt pAt 

= AfiJo dt J M'i^WVt-N 2 , 

(6) 

where ( ) denotes an ensemble average and by the er- 
godic theorem the ensemble average is replaced with the 
time average. Since the correlation function for a Poisson 
distribution is given by (1QJ | 



(N(t)N(t')) t = N 2 + JVe-l*'"*!/ 
the variance becomes 

2Nr 2 



a 2 



At 2 



(e- At / r + At It - 1 



In the limit of At 3> r, we then obtain 

a 2 ~ 2Nt/M 



(7) 



(8) 



(9) 



as expected. 

The probability distribution for Si is then obtained 
from Eqs. (0 and (|Sj) as 



(10) 
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where 




5 = (JVa + b)At 




and the deviation as is giv 


an by 





a s = aaAt -> V2iVsVAt 



(11) 



(12) 



with the arrow indicating approximation under the con- 
dition of At/r — > 00 and T2 = 0. The signal to noise 
ratio (S/N)s for Si is given by 



(S/N)i 



S_ 



(Na + b) 



[N + b/a] 



At 

2Nt 



oc VAi, 



(13) 



increasing as the square root of the bin time. Fig. 
shows the behavior of P(Si) for different ratio of At/r. 
The distribution is single-peaked centered around S and 
the relative width of the peak with respect to the mean 
S becomes narrower as At is made larger, as expected 
from Eq.([T3|). 



B. At < r limit 

In this limit, N(t) remains at a certain integer value 
m through out the bin time and thus 



rU+At ( rn \ 

Si = j t I 5>;(*) + &(*)] 



dt 






where 



■Bi 



U+At 



(14) 



Ajj = / dj(t)dt 

hi 

rti+At 

Bi = b(t)dt 

Jti 



(15) 



are the number of fluorescence and background counts in 
the bin time, respectively, and thus integers. 

In general, the statistics of Ajj and Bi, with associ- 
ated distribution functions P a (Aj yi ) and P^Bi), respec- 
tively, are not necessarily Poissonian. However, in many 
cases these quantities follow Poisson statistics. For exam- 
ple, although the photon statistics of resonance fluores- 
cence of a small number of atoms is sub-Poissonian, when 
measured with an imperfect photodetector, the count- 
ing statistics become Poissonian. In addition, statis- 
tics of scattered light of laser beam is Poissonian. Of 
course, there are cases where these statistics become 
super-Poissonian, particularly when laser power fluctu- 
ations and other technical noises enter. For now, we just 
assume both Ajj and Bi follow Poisson statistics with 
mean values A = a At and B = 6 At, respectively. 

Under this assumption, the conditional probability for 
Si with a constraint N = m is given by 



P(S i \m) = J2 



{Aj,i,B t } 



]\Pa(A^) 



i=i 



Pb{Bi) (16) 



where J2\a ba represents summations to be performed 
for all possible combinations of Ajj and B{ under the 
constraint of Eq. (jTljl . If we assume Poisson distributions 
for P a (Aj,i) and ft (5*), 
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FIG. 2: Signal distribution given by Eq. (fTU)) with iV=1.9, a=3270 s _1 , 6=3140 s _1 and r=23 s. These parameters are derived 
from the experimental data in Ref. [2|. In plotting the distributions, signal counts are divided by corresponding bin times 
for the horizontal axis so as to align the distributions with different bin times. The dotted lines indicate corresponding atom 
numbers. 



where 



mA + B 



(18) 



The resulting distribution is just a Poisson distribution 
with both mean value and variance equal to S m . For 
non-Poisson distributions for P a (Aj^) and P^Bi), the 
resulting P(Si\m) is not Poissonian. However, it is still 
a well-localized Gaussian-like distribution with a mean 
value S m , but its variance is no longer equal to S m . 

The probability distribution for Si for all possible N 
values is then given by 



P(Si)= J2 P ( S i\ N = m ) ' p i 



N= 



(19) 



ra=0 



where each P(S\m), peaked around its mean S m with a 
variance a^, is modulated by P m as shown in Fig. 02 The 
signal to noise ratio (S/N) m for the N = m signal level 
becomes, 



peaks is equal to 



m - 



b/a 



aAt 



(21) 



Unless this ratio is very small, the adjacent peaks sub- 
stantially overlap and thus we have a significant probabil- 
ity of indeterminate atom numbers. A necessary condi- 
tion for well separated adjacent peaks is then r ^> At ^> 
A/a. 

For example, consider the set of parameters used in 
Figs. 2 and 3, a=3270 s _1 , 6=3140 s -1 and r=23 s, but 
with a very short bin time, At/r=0. 00043. For these pa- 
rameters the ratio in Eq. (21) is 0.34, 0.48, 0.59, 0.68 for 
m = 0, 1, 2, 3, respectively, and thereby results in signif- 
icant overlap between adjacent peaks for m > 1. This 
situation is illustrated in Fig. 6(a), which is the result of 
Monte Carlo simulation obtained for this set of param- 
eters. Detailed discussion on Fig. 6 will be given in the 
next section. 



(S/N) r 



(ma + b) At oc V At . (20) 



The half width of the mth peak in the signal distribution 
is also given by the square root of S m and thus the ratio 
of the mth-peak full width to the spacing of two adjacent 



III. NUMERICAL SIMULATIONS 

In the preceding sections, we argued that the two-atom 
loss terms in the master equation are negligibly small 
for a few-atom trap with a few-micron in size and thus 
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FIG. 3: Signal distribution given by Eq. $T$ with JV=1.9, 
a=3270 s _1 , 6=3140 s _1 , r=23 s, At=0.100 s, and thus 
At /r =0.0043. These parameters are derived from the ex- 
perimental data in Ref. 9]. Each peak in the distribution is 
marked with a corresponding atom number m. 



the atom number distribution function is approximately 
Poissonian. When the number of atoms in such micro- 
scopic trap is increased with its size fixed, the two-atom 
loss processes take place more frequently. As a result, 
the atom-number distribution deviates significantly from 
a Poissonian distribution and thus the Poisson approxi- 
mation in the preceding sections are no longer applicable. 
In this section, we include the two-atom loss term and 
calculate distribution functions numerically. 

Although the master equation, Eq. (0), cannot be 
solved analytically, a steady-state solution can be found 
numerically. In steady state, we have dP^ / dt = and 
by rearranging terms we obtain the following recursion 
relation with P_i = 0. 

RP N = ^(N + l)(2T 1 +Nr 2 )P N+1 



1 



(N + 2)(N + l)T 2 P N+2 . 



(22) 



Using this relation the atom-number distribution P/v can 
be easily calculated by iterative method. Alternatively, 
one can calculate a fluctuating time sequence of Si in 
steady state by simulating loading and one- and two- 
atom losses and simulating fluctuating dj(t) and b(t) in 
Monte Carlo simulation. From the time sequence, one 
can calculate the histogram of atom number, i.e., the 
steady-state atom-number distribution. 

We compare the results of these two numerical meth- 
ods in Fig. 21 The values of i?, T\ and Y 2 used in the 
calculations were derived from the experimental data of 
Ref. 9]. A Poisson distribution with the same R and 
Ti is also shown in Fig. 21 (by filled circle-line) for com- 
parison. Once P(N) is known, we can calculate the 
mean atom number N and variance a 2 . The results are 
N = 1.6 and a 2 = 1.5, which should be compared with 
N = a 2 = R/Y 1 = 1.9 obtained for T 2 = 0. With 
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FIG. 4: Atom number distribution for i?=0.080 s~\ Ti=0.043 
s _1 and r2=0.0056 s _1 . For Monte Carlo simulation, we as- 
sumed A£=0.20 s. Empty bar shows the probability calcu- 
lated by iterative method and the filled bar that by Monte 
Carlo simulation. Filled circle-line shows a Poisson distribu- 
tion with the same R and Ti but T2 = whereas empty 
square-line represents a Poisson distribution with the same 
mean value as the iterative and Monte Carlo simulations. 



inclusion of T2 term, the mean atom number decreases 
because of the additional loss term. Although the distri- 
bution is not Poissonian, the deviation from a Poissonian 
distribution with the same TV value is negligibly small. 
According to Ref. 9], the correlation function which in- 
cludes the two atom loss term can be approximated by 
the functional form for Poisson case as in Eq. J7J) with r 
replaced with an effective correlation decay time r e ff . 



Te ff = i/r eff = 4/[Ti + 3^/t 2 + mv 2 ] 

For the above parameters, r e ff=18 s, compared to r = 
17 x =23 s. 

This observation allows us to use Eqs.© and (fTTTj) for 
the calculation of P(Ni) and P(Si), respectively, with 
the N and r e ff values obtained above for nonzero T2. 
For Af < r limit, we can calculate P(Si) distribution 
by using Eq. (|T5)l with substitution of the exact PN=m 
obtained numerically. In Fig. [51 the solid lines are given 
by Eqs. ([TUj) and JTJJ) and the filled area is by the Monte 
Carlo simulation. 

The signal distribution P{Si) in the intermediate re- 
gion, other than two limiting cases considered above, can 
only be obtained by Monte Carlo simulations. From the 
time sequence of Si calculated by means of the Monte 
Carlo simulation with the aforementioned parameters, we 
can calculate P{Si) for various At/r by combining Si val- 
ues in neighboring time bins. The results are summarized 
in Fig. For At/r <C 1, individual atom- number peaks 
are well separated and resolved as shown in Fig. EJb) as 
long as At/r ^> 4/(ar). Otherwise, the peaks for large 
m overlap with neighboring peaks significantly as shown 
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FIG. 5: Signal distribution for i?=0.080 s~\ Ti=0.043 s~\ 
r 2 =0.0056 s _1 , a=3270 s _1 and 6=3140 s" 1 . (a) At > r 
limit, and (b) At <^ r limit. Solid curves in (a) and (b) are 
obtained by Ea. (fTUJ) and Eq.fHIJ), respectively, and filled areas 
are the result by Monte Carlo simulations in both cases. In 
(b), each peak is marked with a corresponding atom number. 



in Fig.EJa). As the ratio At/r increases, the broad back- 
ground appears and grows in height as in Fig. EJc) until 
the background outgrows the atom-number peaks com- 
pletely as in Fig.Efd). 



IV. OPTIMAL BIN TIME 

The trend observed in Fig. HJ can be formulated in a 
quantitative way. We have observed for At C r that 
individual signal distributions significantly overlap with 
neighboring peaks (due to poor signal-to-noise ratio) un- 
less At is much greater than 4/a. The overlap of distri- 
bution functions leads to an increase in the probability 
of having indeterminate atom numbers. We can quan- 
tify this probability f B < as a sum of all areas under the 
distribution function P(Si) outside the boundaries set 
by S m — vA < Si < S m + r]A around the rath peak with 
r] < 0.5. In the time trace picture of fluorescence signal as 
in Fig. 7, this probability is proportional to the number of 
data points outside the region specified by dotted lines 
around a mean signal level. For these data points the 
atom number cannot be assigned unambiguously. From 
Eq. (|19|) we then obtain 



*<(At) = 



1 



1 — 2/7 



i - J2 / p ^ dS * 



(23) 



where the factor 1/(1 — 2t?) is introduced in order to make 
f B < be properly normalized in the limit of 77 — ► 0. 

The atom number also becomes indeterminate if it 
changes during the bin time as in the case of At ^> r. 



From the master equation, Eq. J2J, it can be seen that 
the total rate of change r to t of the atom number is given 

by 



r tot (m) = J R + r 1 m + r 2 



(24) 



for the atom number ra at that instance. The probability 
that the atom number would change from N = ra during 
At is then 



<p{/Ab\N 



ra 



l-exp[-r tot (ra)At]. (25) 



By summing over all possible atom numbers according 
to Pat, we obtain the probability #>(At) that the atom 
number would change during At regardless of its initial 
values. 



*>(At) = J2 i 1 - ex P [-rtotWAt]} • Pi 



N- 



(26) 



A/=0 



If the atom number changes during the bin time At, the 
atom number cannot be determined unambiguously from 
the signal level for this particular bin time. Therefore, 
f P > (Ai) can be regarded as the total probability of inde- 
terminate atom numbers for At > r. 

In general, the above two processes occur indepen- 
dently and thus can occur simultaneous during At. 
Therefore, the total probability of indeterminate atom 
numbers for arbitrary At is given by 

fftot(At) = *>(A*) +2><(At) - *>(At) • <e<{At). (27) 

In Fig. |SJa), this probability $ tot (At) is plotted as the 
ratio At/r for several 77 values. Symbols represent the 
results of Monte-Carlo simulations. The bin time At that 
minimizes this probability can be regarded as an optimal 
bin time for accurate measurement of the instantaneous 
atom numbers in a few-atom trap. The optimal At/r 
value is plotted as a function of 77 in Fig.|HJb). It can be 
seen that for 0.2 < 77 < 0.4 the optimal bin time is within 
the range of 0.003 < At/r < 0.008. Particularly, for 
77 = 0.3, we have an optimal bin time of At/r = 0.004 for 
the same parameter values as used in Figs. 2-6. Among 
the plots in Fig. 6, plot (b) is the most closest to the case 
of the optimal bin time. 



V. COMPARISON WITH EXPERIMENTAL 
DATA 

Detailed information on our experiment can be found 
elsewhere jj, |9(- In short, a few rubidium atoms were 
trapped in a microscopic MOT with a diameter of a few 
microns and fluorescence induced by a trap laser was 
measured in photon counting mode. A raw experimen- 
tal data, a segment of which is shown in Fig. EJa), was 
taken with a bin time of 0.20 s. The atom-number cor- 
relation time r was measured to be 23 s, resulting in 
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FIG. 6: Dependence of signal distributions on At/r. (a), (b), (c) and (d) shows the results of Monte Carlo simulations. 
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FIG. 7: Signal counts outside a region centered around the 
mean signal level cannot be associated unambiguously with 
the atom number corresponding to the mean signal level. 



At/r of 0.0086. Distributions with larger values of At/r 
are derived from the raw data by combining counts in 
neighboring time bins. 

We pay close attention to Fig. 0{b), where the fit is 
given by Eq. (fTU|) with each P(Si\m) given by a Gaussian 
distribution with a mean of S m = mA + B and a variance 
of a^ to be evaluated below. 

The average background photon number, B = b At, 
and the average level spacing, A = aAt, are 627 and 653, 



respectively, obtained from the experiment. By counting 
the individual loading and loss events in the time trace 
of fluorescence as shown in Fig. 0{a), one can measure 
the loading rate R and the one- and two-atom loss rates 
Ti and T2, respectively, and the results are i?=0.080 s _1 , 
Ti=0.043 s" 1 , r 2 =0.0056 s" 1 . The detailed information 
on experiments to measure these rates can be found else- 
where 9]. From these parameters we obtain N= 1.6. The 
Poisson distribution for this actual N is used for P/v= m 
in Eq. (|19|) for the fit. Note that the only fitting param- 
eter is then the variance cr^, which can be decomposed 
into 



2 2 

^m =(7 B 






(28) 



where <j b and <j| are the variances of background and 
signal counts, respectively. 

The background counts are mostly due to scattered 
light of trap and repump lasers of MOT. Due to long- 
term power fluctuations, the mean value of background 
counts also fluctuates, and as a result, the width of the 
zero- atom peak in the signal distribution P(Si) becomes 
larger than that of a Poissonian distribution. In fact, the 
background variance <j b was measured to be 2570±180, 
about 4.5 times larger than the mean count B. 

If we assume that the fluorescence counts follow Pois- 
son statistics, the variance <r^ can be modeled as 



0"™ = 



4 



ma \ = (2570 ± 180) + mA (29) 
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FIG. 8: (a) Total probability ^Ptot(At) of indeterminate atom 
numbers, given by Eq. (|27f> . for various values of r\. Symbols 
denote Monte-Carlo simulation results, (b) Optimal bin time 
as a function of r\. 



where a\ is the variances of one- atom fluorescence, and 
it is assume that the fluorescence from one atom is statis- 
tically independent from that of another atom. However, 
as shown in Fig. Ofc), the observed variances of individ- 
ual peaks are not well fit by the above formula. Rather 
they are well fit by an empirical formula given by 



a^ = (2570 ± 180) + (2280 ± 190)ra 



(30) 



the slope of which is about four times larger than that of 
Eq. (2J . 

The fact that the variance is still linear in m indicates 
that the fluorescence from one atom is still statistically 
independent from that of another atom. This observation 
excludes, as a source of the increased variance, the fluo- 
rescence dispersion due to power fluctuation of trap and 
probe lasers, mechanical vibrations and similar technical 
noises since they all have to induce correlated fluctuations 
in the signals of individual atoms and thus proportional 
to m 2 . 

One possible reason for this increase variance is the 
motional effect of individual atoms. The atoms move 
independently from each other inside the MOT. Due to 
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FIG. 9: (a) A segment of the fluorescence signal time trace 
observed in the experiment in Ref. |9|. Loading events and 
atom- and two-atom loss events can be identified in the trace, 
(b) Signal distribution obtained from the experimental data in 
(a). The filled area shows the experimental result. The solid 
line shows a fit based on Eqs. (fT§J) and (|30J). Each peak is 
marked with a corresponding atom number, (c) The filled 
squares represent the variances of individual atom-number 
peaks in the observed signal distribution in (b) whereas the 
unfilled squares show those of Poisson statistics. Bin time At 
was 0.20 s and the correlation decay time r was about 23 s 
in both experiment and analysis. All other parameter are the 
same as in Figs. 2-6. 



the spatially inhomogeneous magnetic field, atoms expe- 
rience different Zeeman shifts and thus their upper level 
populations vary in time differently and independently 
from one atom to another. This variation can give rise 
to the observed increased variance in fluorescence counts. 
Such motional effect might be observed in the second 



9 



order correlation function of the fluorescence in the long 
time limit. In the short time limit, comparable to the 
life time of the atom (tens of nanosecond) , antibunching 
characteristics of the resonance fluorescence will be dom- 
inant effect. But in the long time limit, much longer than 
the atomic life time and comparable to the characteristic 
time (~ millisecond) of atomic motion in the trap, an 
oscillatory feature would appear in the second order cor- 
relation function. The detailed study on this phenomena 
is beyond the scope of this paper and left for the future 
work. 



VI. CONCLUSION 

We have derived analytic expressions for signal distri- 
bution P(Si) of fluorescence photo-counts from a few- 
atom MOT and compared the results with Monte Carlo 
simulations and experimental data. The signal distri- 
bution strongly depends on the relative size of the bin 
time At of photon counting with respect to the trap de- 
cay time r. In the limit of At <C r, the distribution 
shows multiple peaks with the integrated areas of indi- 



vidual peaks constituting the atom-number distribution 
function P/v- Conversely, the stepwise fluorescence sig- 
nal corresponding to a multi-peak distribution can be 
regarded as a definitive evidence of a few atoms in the 
trap. As At C r is increased, a broad background ap- 
pears and eventually outgrows sharp peaks corresponding 
to atom numbers and turns into a single peak in the limit 
of At ^> r. The validity of our derivation was confirmed 
by comparing the results with those of numerical simula- 
tions including Monte Carlo simulation. These theoreti- 
cal results were then compared with experimental results. 
Fluorescence photo-count distributions were observed to 
be super-Poissonian, the origin of which might be due 
to the statistically independent motion of atoms in the 
inhomogeneous magnetic field of MOT. Our results pro- 
vide necessary theoretical background for analyzing and 
interpreting the fluorescence signal of a few atom MOT 
and also clarify the optimum condition on the bin time 
in actual experiments. 
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